DFT studies of ethylene in femtosecond laser pulses 
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Abstract 

Using time-dependent density functional theory, apphed to valence electrons, coupled non- 
adiabatically to molecular dynamics of the ions, we study the induced dynamics of ethylene sub- 
jected to the laser field. We demonstrate the reliable quality of such an approach in comparison 
to the experimental data on atomic and molecular properties. The impact of ionic motion on 
the ionization is discussed showing the importance of dealing with electronic and ionic degrees of 
freedom simultaneously. We explore the various excitation scenarios of ethylene as a function of 
the laser parameters. We find that the Coulomb fragmentation depends sensitively on the laser 
frequency. The high laser intensity can cause brute-force Coulomb explosion and the laser pulse 
length actually has influence on the excitation dynamics of ethylene. 
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I. INTRODUCTION 



The developments in laser technology have opened the door to many new research areas 
of physics, such as molecular physics. It is a powerful tool to trigger and to analyze a 
bunch of dynamical processes of molecules and clusters. These processes include various 
dynamical situations, such as optical response [l], multi-photon ionization j^, and Coulomb 
explosion of clusters . When the intensity of laser is comparable to molecular electronic 
fields, the violent excitations quickly drive the system beyond the Born-Oppenheimer surface 
resulting some new processes such as high-order harmonic generation 4,l5|,l6|, bond softening 
0, [sl, charge resonance enhanced ionization [o], 1^. The theoretical understanding of these 
mechanisms needs to solute the time-dependent Schrodinger equation for all electrons and 
all nuclear degrees of freedom. However, this kind of numerical solutions only exist for 

12], H2 and H+ Q, 



smallest systems, such as atoms 



11 



14 



15 



llTj . For large systems 



the computation is quite expensive due to the number of degrees of freedom. 

Recently, Horsfield et al. use correlated electron-ion dynamics to study the excitation of 
atomic motion by energetic electrons [l^. Saalmann et al. have developed a non-adiabatic 
quantum molecular dynamics (NA-QMD) to study different non-adiabatic processes of differ- 



ent systems 
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of rare-gas clusters 



201 . 21 1 and Calvo et al. use a combined methods to study the fragmentation 



22l |. In contrast, a fully-fiedged coupled ionic and electronic dynamics is 



developed. At the side of the electronic dynamics, the well tested time dependent Density 
unctional Theory at the level of the time-dependent local-density approximation (TDLDA) 



23( 1 is applied. The ions are treated by classical molecular dynamics. This altogether pro- 



vides a coupled TDLDA-MD, for a review see 



241 ]. Up to now, it has been applied to study 
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3l|. 



the dynamical scenarios of simple metal and hydrogen clusters 
Besides these, organic molecules are particularly interesting cases now, such as ethylene. 
For ethylene, recently a lot of theoretical work has been carried out to investigate either the 



the external laser field 
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optical response in linear domain 32| or the dynamics of it without and with considering 



33|. 



In this paper, applying TDLDA-MD, we first devote to exploring the impact of ionic 
degrees of freedom on the excitation dynamics and we find the enhancement of the ionization 
with ionic motion. Furthermore, we study the various excitation scenarios of ethylene as a 
function of laser parameters. We find that these scenarios depends sensitively on the laser 
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parameters. It should also be noted that switching from simple metals and hydrogen clusters 
to organic molecules is a challenging test for the model. 

The paper is outlined as follows. Section 2 provides a short presentation of the theo- 
retical and numerical approach. Section 3 first gives the ground-state properties of atoms 
and molecules. Then the impact of ionic motion on the excitation dynamics is discussed. 
Finally, the different excitation scenarios of ethylene as a function of the laser parameters 
are presented. Section 4 gives conclusions. 



II. THEORETICAL AND NUMERICAL METHODS 



The technicalities of the model we use have been presented in detail elsewhere and for a 



review see 



24l ]. Here for the sake of completeness, we recall in this section the ingredients 



and a few relevant formulas. 

The molecule is described as a system of valence electrons and ions. Valence electrons 
are treated by TDLDA, supplemented with an average- density self-interaction correction 
(ADSIC) 3^. Ions are treated as classical particles. They interact via their repulsive 
Coulomb force. Pseudopotentials are used for the interaction between ions and electrons. 
We apply the form from 35|. The TDLDA equations are solved numerically on a grid 3D 
coordinate space (cuboid box size of 72*72*64) with a grid spacing of A = 0.41 ao. 

The computation of ground state wavefunctions is determined by the damped gradient 
method 36|. The time-splitting method is applied to solve the time-dependent Khon-Sham 



equation for electrons. For ions a Verlet algorithm is used for their propagation under their 
mutual Coulomb interactions together with the electron-ion forces as driven from the pseudo 
potential. For the nonlocal part in the pseudopotentials, it is treated as a third-order Taylor 
expansion of the exponential 37|. The absorbing boundary condition is employed to avoid 
periodic re- feeding of the emitted electrons 38[ . For the laser field, neglecting the magnetic 
field component, acting both on electrons and ions, can be written as 



V, 



laser 



EoZf laser (t) COS (uJlasert) 



(1) 



where Eq oc vT, / denoting the intensity, z is the dipole operator, fiaser{t) is the pulse profile 
and ujiaser IS the laser frequency. In this paper, fiaserif) is chosen as cosine^ in time. It should 
be noted that this coupled ionic and electronic dynamics, constitutes a true TDLDA-MD, 



3 



and goes beyond the usual Born-Oppenheimer MD. It is powerful to deal with the violent 



'he capability of this technique 
and hydrogen clusters [30 1. 



off-equilibrium dynamics as induced by strong laser pulses, 
was demonstrate earlier calculations both on metal clusters 
For organic molecules, in order to get the confidence of our approach to describe them, we 
first start with the ground state calculation of ethylene. Here we devote to the ionization 
potential (IP) of Hydrogen, Carbon and ethylene. The IP which is important for dynamical 
calculations is computed from the energy of the last occupied electron state. In addition 
we calculate the bond lengths of ethylene. For the more dynamical observables we study 
excitation spectra of ethylene in the frequency domain which is quite useful for exploring 
the excitation under the laser pulse. The excitation spectra we get is quite similar to what 



Takashi et al 



32l | got. To avoid the repetition, in this paper we give the location of the 



lowest peak and compare it with the experimental data. 



III. RESULTS AND DISCUSSION 



A. Ground state properties 




FIG. 1: Ionic structure of ethylene 



The ethylene, C2H4, is the simplest organic vr system holding symmetry. There are 
12 valence electrons. In our calculation, it is placed in x-y plane and the center of mass is 
at the origin. The ionic structure of C2H4 is displayed in Fig. 1. 

Table 1 shows the ionization potentials (IP) of neutral atoms and C2H4. It is obvious 
that the IP are reproduced well when compared with experimental values. Table 2 gives 
more detailed results of ethylene. One can find that the relative errors are less than 4% and 
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TABLE I: Ionization potentials (IP) for H, C and C2H4 obtained with LDA-SIC. In unit of eV. 



Experimental data are from the website 



39|. 



H C C2H4 



exp [eV] 13.6 11.3 11.0 
LDA-SIC [eV] 14.0 11.3 11.5 
% error 2.9 4.5 



TABLE II: Theoretical bond lengths, in unit of a.u., and electronic excitation energies of C2II4, in 
unit of eV, acquired from LDA-SIC. Experimental data come from [39, 40]. 





Rcc[a.u.] 


RcH[a.n.] 


Lowest Peak [eV] 


exp 


2.53 


2.05 


7.66 


LDA-SIC 


2.43 


1.97 


7.74 


% error 


4.0 


1.5 


1.0 



the theoretical results are quite acceptable. 
B. Dynamics driven by laser pulse 

The laser field is coupled to ethylene as a time-dependent classical field interaction simply 



via the dipole operator. As is discussed in 4l|, the laser pulse is characterized by frequency, 



intensity and width which have their own specific infiuence. The frequency determines 
the strength of the laser-molecule coupling and the intensity determines the dynamical 
regime. The pulse length determines the spectral resolution and the interplay with ionic 
motion. Thus it is painful to choose the test case in the enormous variety of conceivable 
laser parameters. It should be noted that in all cases discussed following the laser pulses are 
switched off at the moment of t=Tpuise- 

1. impact of ionic motion 

To investigate the effect of ionic motion on the excitation dynamics, we compare cal- 
culations either with full ionic motion or with fixed ionic configuration. Fig. 2(a) and (b) 
show the time evolution of ionic motion and the number of escaped electrons (Ne^c) with full 
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FIG. 2: A comparison of the excitation of ethylene induced by a laser pulse with full ionic motion 
(a) and with fixed ionic configuration (b). From the top to the bottom panels in both parts are time 
evolution of ionic position along laser polarization and the number of escaped electrons (carbon: 
hollow square, hydrogen: full square). The laser polarization is along x axis. I = 10^^ W/cm^, 
i^iaser = 8.16 cV and rp„/se=100 fs. 

ionic and with fixed ionic motion respectively. The difference between two cases is dramatic. 
From Fig. 2 (a) we can see that ions start to move at around 40 fs. Four hydrogen ions move 
more obviously than carbon ions. Electrons start to escape gradually at around 18 fs. At 
around 50 fs, when the laser reaches its highest amplitude, the electron emission increases 
rapidly and there are almost 2.8 electrons emitted when the laser is switched off. In Fig. 2(b), 
when ions are kept fixed during the calculation, the trend of electronic emission in the first 
40 fs is almost the same as in Fig.2(a). However, the difference of Nesc between two cases 
appears at around 40 fs, which is the same moment at which ions start to move. Finally 
when the laser is switched off, about two electrons are emitted. This value is smaller than 
that of the moving ions case. The enhancement of the ionization for the moving ions case 
can be explained as in [^^J even if the studied system here is the covalent molecule, ooiaser 
is in the resonant region, after increasing ionization, the optical spectrum has blue shift, 
so uiaser will be Smaller than the resonant frequency. However, the Coulomb force from 
increased net charge impels the expansion of ions which could make red shift of the optical 
spectrum and Uiaser becomes in the region of resonance again. 

The difference of electron emission above shows the importance of treating electronic and 
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ionic dynamics simultaneously. With the coupling of electronic and ionic degrees of freedom 
we can explore whatever excitation regimes. 

2. influence of laser frequency 




20 40 60 80 20 40 60 80 
Time [fs] Time [fs] 



FIG. 3: Excitation of ethylene by two types of laser pulses with different frequences, (a) ujiaser = 
9.5 eV, (b) uJiaser = 6.8 eV. / = 10^^ W/cm^ and Tp^ise = 20 fs. The polarizations of lasers are 
both along the x axis. From the top to the bottom panels in both parts are the time evolution of 
bond lengths, extensions of the ionic distribution in three directions, dipole signal along the laser 
polarization and the number of escaped electrons. 

We now turn to the influence of the laser frequency on the excitation dynamics of ethylene. 
We consider two frequency cases here, one is uiaseT = 9.5 eV (Fig. 3(a)) which is in the 
resonant region and the other one is ujiaser = 6.8 eV (Fig. 3(b)) which is off the resonant 
region. From the top to bottom panels are the time evolution of the bond lengths, of the 
ionic extensions along x, y, z directions, which are ^/{x^, (y^), ^y {z^), of the dipole signal 
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along the laser polarization, and of l^esc- The polarizations of laser pulses are both along x 
axis and / = lO^'^ W/cm^. As the Tpuise = 20 fs is rather short, after laser is switched off, 
we keep on following the tracks of ions and electrons until 100 fs. 

From the top panel in Fig. 3(a), one can see that all of CH bond lengths start to increase 
slowly in parallel at around 10 fs. After the laser is switched off, they keep on increasing 
and start to separate from each other at around 25 fs. All of them are larger than 10 a.u. 
at around 60 fs. Such a large bond length reflects that CH bonds are broken. The CC bond 
length (solid line) starts to change a little later than that of Kch and oscillates gently during 
the whole time without broken. For the ionic expansion, as is shown in the second panel in 
Fig. 3(a), it is obvious that ethylene expands within its plane and there is no expansion in z 
direction (dotted line). Compare the top two panels in Fig. 3(b) with corresponding panels 
in Fig. 3(a), we can find that for the off-resonant case, ions stay well around their initial 
positions. 

More physical insight can be gained by looking at the dipole moment as shown in the 
third panel in both parts. When ujiaser = 9.5 eV, it has no resemblance with the laser pulse 
profile. The first obvious burst appears at around 10 fs. This is almost the same moment 
at which the bond lengths start to increase and the laser reaches its highest amplitude. At 
around 20 fs, when the laser is switched off, the dipole signal arrives its highest amplitude 
and there are many bursts. It keeps on oscillating strongly after the laser is switched off. 
At around 58 fs, a quick shift of the dipole signal takes place. This is due to the fact that 
hydrogen ions are moving towards the border of the calculation box. In contrast, for the 
off-resonant case, the dipole signal has almost the same profile as the laser pulse and it 
attains highest amplitude when the laser field reaches its highest field. After the laser is 
switched off it remains oscillating gently. 

The dipole amplitude has consequence for the direct electron emission, which is the same 
as in the cases 



strong emission while small amplitudes cause weak emission. As we can see the electron 
emission in the bottom panels. In the bottom panel in Fig. 3(a), the N^sc starts to increase 
at around 5 fs with a constant speed. This moment is a little earlier than the movement of 
ions since electrons are more active than ions. When the laser is switched off at 20 fs, there 
are almost 2.5 electrons escaped. The ionization continues but becomes much more slowly 
after the laser is switched off. However, at around 58 fs, ionization increases quickly again. 
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For the electron emission in the bottom panel in Fig. 3(b), it starts to increase at around 
5 fs and levels off quickly after the laser is switched off at 20 fs. The final value of Nesc is 
much smaller than that in Fig. 3(a). 

From above discussion we can recognize that the laser frequency plays an important role 
in the excitation dynamics of ethylene. The frequency in the resonant region enhances the 
ionization leading to the Coulomb fragmentation. 

3. influence of laser intensity 



l=10''W/cm^ 
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FIG. 4: The time evolution of ionic motion along the x direction and ionization of ethylene with 
different laser intensities. The laser polarization is along the x direction, Tp^igf. = 40 fs and 

^laser = 6.8 eV. 

We now keep the laser frequency and the pulse length fixed and vary the laser intensity. 
Fig. 4 shows the excitation dynamics of ethylene by the laser pulses with different intensities. 
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The ujiaser IS 6.8 eV which is off-resonant and Tpuise = 40 fs. One can find that when 
intensities are low (left panels), ions are almost frozen in the course of the laser pulses and 
electron emissions in the bottom panel are quite small. Ne^c values are multiplied by 1000 
and 50 respectively in the bottom panel in the left part in order to compare with Ng^c 
values in the bottom panel in the right part. In the right part, for I = 10^^ W/cm^, ions 
seem to move a little at the end of the laser pulse. The dissociation of ethylene occurs 
when / = 10^^ W/cm^. Ions start to move at around 10 fs with a constant velocity. It is 
clear that hydrogen ions move faster than carbon ions. From the right bottom panel we 
can see that when I = 10^^ W/cm^, eight to nine out of the twelve electrons are emitted 
in a comparatively short time. The ethylene is grossly Coulomb unstable at this charge 
state. As a consequence, the ionic motion exhibits clearly the pattern of a straight Coulomb 
explosion with CC double bond and CH bonds broken as shown in the middle panel of the 
right part. Comparing the electron emissions of different intensities, one can also find that 
for I = 10^^ W/cm^, the Nesc (dash-dotted line) starts to increase at around 7.5 fs. This 
moment is earlier than not only the ionic motion but also that of other cases with lower 
laser intensities. 

The results above show that the ionization of ethylene is enhanced by the high laser 
intensity and the high intensity can cause Coulomb explosion even if the laser frequency is 
in the off-resonant region. Ions exhibit rigidity when laser intensities are low. 



4- influence of pulse length of laser 
Christov et al 



431 ] found that the harmonic spectrum of argon taken with 25 fs laser 
pulses contains harmonics up to 20 orders higher than for 100 fs laser pulses with the same 
intensity due in part to the nonadiabatic response of the atomic dipole to the fast rise time 
of the pulse. Here we explore the influence of pulse lengths on the excitation of ethylene. 
Results are shown in Fig. 5. We present two cases with different laser frequencies, one is 
in-resonant, oJiaser = 13.6 eV (Fig. 5(a)) and the other one is off-resonant, ujiaser = 6.8 eV 
(Fig. 5(b)). We choose the moderate peak intensity in both cases, / = lO'^^ W/cm^. Since 
durations are short for Tp^ise = 20 fs and 40 fs, after the laser is switched off, we keep on 
tracking the dynamics of electrons and ions until 60 fs and 80 fs, respectively. 

In Fig. 5(a) one can find that with the same intensity, when frequency is in the resonant 
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FIG. 5: Influence of laser pulse lengths on the excitation of ethylene. Tp^^gg = 20 fs, 40 fs, 100 fs. 
The laser polarizations are along x direction and / = 10^^ W/cm^. The time evolution of ionic 
motion along the laser polarization and ionization are presented. Two types of laser frequencies 
are given, (a): ijJiaser = 13.6 eV, (b): ujiaser = 6.8 eV. 



region, no matter how long the pulse length is, ions move visibly and hydrogen ions move 
stronger than carbon ions. The difference lies in that for short pulse lengths, ions move 
gently in the course of laser pulses and the post-laser ionic motion are also gentle. However, 
when Tpuise = 100 fs, after the first gentle oscillation, ethylene is dissociated with CH bonds 
broken before the pulse is over. Therefore in this case, the pulse length actually has strong 
influence on the ionic motion and the long pulse length can cause fragmentation of ethylene. 
For the electron emission, one can read off in the bottom panel that electrons start to escape 
earlier than ionic motion and the ionization is enhanced by the long pulse length. This is 
due to the fact that with the long pulse length, active electrons absorb energy from the 
laser field and multielectron excitation leads to a highly excited ethylene with subsequently 
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fragments, assisted by Coulomb repulsion. 

When uiaser is in the off-resonant region, as shown in Fig. 5(b), the pulse length has 
little effect on the ionic motion. In top three panels, ions are rigid in the courses of the 
laser pulses. The time evolutions of Ne^c are shown in the bottom panel. In order to see 
curves clearly we multiply Ng^c values by 20. It is clear that the ionization is enhanced by 
the long pulse length. However, for different laser pulse lengths, ionizations are quite small. 
Furthermore, one can find that the electron emissions are saturated before the laser pulses 
are switched off, no matter how long the pulse length is. The independence of the excitation 
of ethylene on the laser pulse length here is related to the fact that the laser frequency is 
off-resonant and laser intensity is moderate, as a consequence, the laser pulse is not strong 
enough to remove electrons and there is no obvious ionic motion. 

We can learn from the above discussion that when the laser pulse is with moderate 
intensity and in-resonant frequency, the excitation dynamics of ethylene depends sensitively 
on the pulse length. 

IV. CONCLUSIONS 

In this paper we have demonstrated the capability of TDLDA-MD to describe the molecu- 
lar dynamics and the ionic motion really has strong effect on excitation dynamics. We have 
shown that this approach produces reliable basic atomic and molecular properties when 
compared with experimental values. We discussed various excitation scenarios of ethylene 
subjected to different laser pulses. These scenarios involve both electrons and ions but the 
relative role of each species does depend on the actual excitation conditions of the laser 
pulse. First, varying the laser frequency with fixed laser intensity and pulse length indicates 
that the appearance of deexcitation/explosion scenarios depends on the relation of the laser 
frequency to the eigenfrequencies of the system. Second, with fixed laser frequency and pulse 
length, changing the laser intensity reveals that the high laser intensity leads into the the 
brute-force Coulomb explosion of ethylene. Finally, we find that with the same moderate 
intensity, the time scales of electronic and ionic motions depend strongly on the pulse length 
when laser frequency is in the resonant region. In contrast, the excitation of ethylene is 
much less effected by the pulse length when the laser frequency is off-resonant. 
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